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Cd Abstract 



Linear singularly perturbed convection-diffusion problems with characteristic 
layers are considered in three dimensions. We demonstrate the sharpness of our 
recently obtained upper bounds for the associated Green's function and its deriva- 
tives in the Li norm. For this, in this paper we establish the corresponding lower 
bounds. Both upper and lower bounds explicitly show any dependence on the sin- 
^ gular perturbation parameter. 
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1 Introduction 

o 

CN Consider the convection-diffusion problem in the domain = (0, 1) : 

in 

Cxu{x) = —e/S.xu{x) — 2adx^u{x) = f{x) for a:; G 17, (la) 

g u{x) =0 for ic G dVL. (lb) 

Here e E (0, 1] is a small positive parameter, while a is a positive constant. Then ([T]) 
ju^ is a singularly perturbed convection-dominated problem, whose solutions typically 

exhibit sharp characteristic boundary and interior layers. 
^ This article addresses the sharpness of our recently published obtained upper 

d bounds for the associated Green's function and its derivatives in the Li norm. 

Our interest in considering the Green's function of problem is motivated by the 

numerical analysis of this computationally challenging problem. More specifically, 
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Figure 1: Anisotropy of the Green's function G{x; ^) associated with ([T]) for e = 0.01 
and X = (|, ^, |). Left: isosurfaces at values of 1, 4, 8, 16, 32, 64, 128, and 256. Right: a 
two-dimensional graph for fixed ^3 = X3. 



these estimates will be used in the forthcoming paper [3] to derive robust a posteriori 
error bounds for computed solutions of this problem using finite-difference methods. 
(This approach is related to recent articles [101 E], which address the numerical 
solution of singularly perturbed equations of reaction-diffusion type.) In a more 
general numerical- analysis context, we note that sharp estimates for continuous 
Green's functions (or their generalised versions) frequently play a crucial role in a 
priori and a posteriori error analyses [3l [U [11] . 

For each fixed x & Q, the Green's function G associated with ([T]) satisfies 

£lG{x;^) ■.= -eA^G{x;i,) + 2adi:,G{x-$,) = S{x-$,) for $ en, (2a) 

G{x;$,) = ioT$£dQ. (2b) 

Here is the adjoint differential operator to Cx, and 5{-) is the three-dimensional 
Dirac (5-distribution. 

Note that the Green's function for a singularly perturbed self-adjoint reaction- 
diffusion operator —eAx -|- a is almost radially symmetric and exponentially de- 
caying away from the singular point [2]. By contrast, the Green's function for our 
convection-diffusion problem ([T]) exhibits a strong anisotropic structure, which is 
demonstrated by Figure [T] 

In 1^1 [n] we have obtained certain upper bounds for the Green's function associ- 
ated with a variable-coefficient version of which we now cite. 

Theorem 1 (|S1[S])- Let e G (0,1]. The Green's function G associated with ([T]) on 
the unit cube Q = (0, 1)^ satisfies, for all a; E $7, the following upper bounds: 

||%G(a;;-)||i;Q<C(l + |lne|), (3a) 

||%G(a;;-)||i;Q < k = 2,3. (3b) 

Furthermore, for any ball B{x',p) of radius p centred at any x' £ Q, we have 

< Cp/e, (3c) 
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while for any ball B(x,p) of radius p centred at x £ we have 



\\dlG{x;-)\\i.n\B(..,p) < C^-^ ln(2 + e/p), (3d) 
\\dlG{x;-)\\^.u\B(.^p)<Ce~Wlne\+ln{2 + e/p)), k = 2,3. (3e) 

Remark 2. Theorem [T] is given in [5, 6j for a more general variable-coefficient op- 
erator —eAx — a{x) + b{x) with sufficiently smooth coefficients that satisfy 
a{x) > 2a > and b{x) + dxia(x) > for all x £ Cl. 

The purpose of this paper is to show the sharpness of the bounds of Theorem [T] 
up to a constant e-independent multiplier in the following sense. 

Theorem 3. Let e G (0,co] for some sufficiently small positive cq. The Green's 
function G associated with the constant- coefficient problem ([T]) in the unit cube 
n = (0, 1)^ satisfies, for all x £ |]^, the following lower bounds : 

\\d^,G{x;-)\\i;n>c\lne\, (4a) 
\\d^^G{x;-)\U.n>ce-^/^ k = 2,3. (4b) 

Furthermore, for any ball B{x; p) of radius p < I, we have 



cp/e, ifp<2e, 
c^p/e)^^"^, otherwise, 
\\dlGix;-)\\i.n\B{a^,p)>ce-Hn{2 + £/p), ifp<ci£, (4d) 

\\dlG{x;-)\\^.n\B(.,p)>ce-\ln{2 + e/p) + \lne\), ifp<\, (4e) 

where k = 2,3, and ci is a sufficiently small positive constant. 

Remark 4. The restriction x £ [j; f]'^ is somewhat arbitrary and can be replaced by 
X G [9,1 — 9]^ with any e-independent constant 9 G (0, ^) (then p < ^9 is imposed 
instead of p < g). 

The paper is structured as follows. Sharp lower bounds for the fundamental 
solution in M'^ are derived in Section [2j The next Section [s] is devoted to the proof 
of Theorem [3] in a simpler setting, for the domain (0,1) x M^. Sections |4] and [5] 
briefly describe an extension of the analysis of Section [2] to the domain Q = (0, 1) 
and to convection-reaction-diffusion problems. Finally, in Section [6j we give an 
outlook for problems in n dimensions. 

Notation. Throughout the paper, C denotes a generic positive constant, typically 
sufficiently large, while c denotes a sufficiently small generic positive constant; they 
take different values in different formulas, but are independent of the singular pertur- 
bation parameter e. The usual Sobolev spaces W"^'P{D) and Lp[D) on any measur- 
able set D C are used; the Lp{D) norm is denoted by ||-||p;D) while the W^'^{D) 
norm is denoted by ||-||m,p;D- By x = {xi,X2, X3) we denote an element of M'^. For an 
open ball in M?, we employ the notation B{x' , p) = {x ^M? : Yl\=i{^k—x'^Y < P^}- 
The notation dxif, d^^f and Ax is used for the first- and second-order partial deriva- 
tives of a function / in variable xi, and the Laplacian in variable x, respectively. 
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2 The fundamental solution 



In this section we investigate the fundamental solution g that solves a similar prob- 
lem to ([2]) but posed in the domain M?: 

Clg{x-$,) = -e^^g{x-^) + 2ad^,g{x-$,) = 5{x-i) for ^ G Ml (5) 

Using [2j ,9j (see also [5j,t6j), a calculation shows that the fundamental solution g is 
explicitly represented by 



9 = 9{x]C) 
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(6) 



with the scaled variables ii^[x^] = 6 = ^3 = and the scaled 

distance between x and ^ denoted by ry^^] = [^i] + ^ + ^he subindex 

[xi] in ^iji^i] and f^xi] to highlight their dependence on xi as in many places xi will 
take different values; but when there is no ambiguity, we shall simply write ^1 and 
r. 

Next, we evaluate derivatives of g of order one and two: 



r 



ga(gi-r) 

ga(fi-r) 
ga(fl-f) ^ 

r 



%5 

dig - 



a 



(af + 1) l|. A; = 2, 3, 



a r 



A; = 2, 3. 



(7a) 

(7b) 
(7c) 

(7d) 



Lemma 5. Let 0* := (0, 1) x M? , e G (0, cq] /or some sufficiently small constant cq, 
and < a < C. T/ien i/ie function 5 0/ Q satisfies, for any x € [2e, |] x M^, ^/le 
following bounds 



\\dii9{x;-)\\i;n, > c\lne\, 
||%5(^c;-)lli;f^, >ce-i/l A: = 2, 3. 

Furthermore, for any ball B{x; p) of radius p < I, we have 



(8aj 
(8b) 



II ^ Ml ^ jcp/e, if p < 2e, 

\\9[x;-)\\i,i-n,nB(x-p)^\ , , si/2 ^i, 

I c(p/e) ' , otherwise, 

\\dlgix;-)\\i.Q^\Bix;p) > ce"Mn(2 + e//>), if p < cie, 

\\dl9{x;-)\\i.nAB(x;p)>ce-HH'^ + e/p) + \lne\), ifp<l, 

where k = 2,3, and c\ is a sufficiently small positive constant. 



(8c) 

(8d) 
(8e) 
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Remark 6. The statement of the above Lemma [5] is almost identical with the state- 
ment of Theorem |3] with G replaced by g and ri* replaced by fi. 

Proof. Throughout the proof, xi is fixed, so we employ the notation := ^ij^i] 
and ^ := (^1,^2,^3)- First, we rewrite all integrals that appear in ^ in variable ^ E 
0^, := (— ^-^^) X using = e^d$. Now it suffices to prove the desired lower 
bounds on any sub-domain of 17*. In particular, we employ the non-overlapping 
sub-domains Jli and 0,2 of fi*: 



Qi := B{0, 1) n {e? > i + el}, ^2 := {max{l, ^^1 + ^1 < 6 < 

(but other sub-domains of B{0, 1) similar to will be considered as well). The 
notation [v]~^ := max{?;,0} will be used for any function v. 

(i) The bounds (8a), (8b) will be obtained using il.2- Note that for $, £ ^2 one 

has <r < -v/2 Ci- Introduce the new variables := ^k/y 2^1 for k = 2,3, so 



where C2 := 2(1 + \/2)~^, and we used f — = (^ + ^)/('^i + r) to get the final 
relation above. 

Now a calculation using (7a) yields the first desired bound (8a) as follows: 

6 



||%5lli,n*> c 




>c er^d^i 

> c I In el. 



air - ^ij 



A similar calculation using (7b) with k = 2,3 yields (8b); indeed, 

3"(6-)(af+l)er='|e,|d^ 



|%<?||l,f2.> C 



e 



(ii) To show (8c) for p < 2e, we note that ||fi'||i,i;B(a;,p) > l|f^65lli;B(a;,p) ^° 
p := pje and consider the sub-domain $73 :=i?(0,p)n|^ >Ci+l3|- Note that 



in this sub-domain, e"*-^^ > c and ^2 ^ ^I\f2 so (7b) yields 

ll%fi'lli,n.nB(£c,p)> c /"/"/ f^^d^>c / dr>cp = cp/e. 



which immediately implies (8c) for p < 2e. 
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Next, for p G [2e, |] consider d^^g in the sub-domain 1^2 H 5(0, p). Imitating the 
calculation in part (i), one gets 



l%5lll,n.nB(a;,p) > C 




■-3/2 



e 



'Q2nB(o,p) 
> c ^dCi (a^i + 1) Ci 

which completes the proof of ( 8c ) for p < | . 

(iii) To obtain (8d), we use the sub-domain Oi\S(0,p), where p := pje. Note 
that for ^ G r^i one has e^'^^^''-' > c and -^^2 — > 1. So using (|7c|), one gets 




> ce 



-1 



0-4af+l di 

Qi\B(0,p) 
min{l, 

r"^ df > ce"-^ ln(2 + p/e). 



> 



So we have shown (8d) for p < cie with ci := | min{l, g^}. 

(iv) In a similar manner as in part (iii), using (7d), one can show that ||9|^5'||i,Q\B(a;,p) 
ce~^ln(2 + p/e) for A; = 2,3 and p < ^e. Note that now we use the sub-domain 
B{0, 1) n > ^ + ^1 instead of hi, with j = 3 for A: = 2 and j = 2 for A; = 3. 

Consequently, to obtain (8e) for any p < |, it remains to show that ||c^|^5||i,r2*\_B(a3,p) ^ 
ce~^\ lne|. For this, consider (7d) in il,2\B{0,p). Combining the observations made 
in part (i) with ar + 1 < (a\/2 + and 



> — 1, yields 



\\dl^9\\i,n.\B{x,p) 



> ce 



> c 



-1 



^«(^i-f) ?-3 



max{l,p} 

> ce~^| Inel. 



n2\B(0,p) 



2aVfc - (a\/2 + 1] 



Here we also used p < |e So we have proved the final desired assertion (8e). □ 

□ 



3 Approximation of the Green's function and 



proof of Theorem 3 for the domain = (0, 1) x 

To approximate the Green's function, we use the fundamental solution g of Section[2] 
and the cut-off function w, defined by 

uj{t) e C^{0, 1), u>{t) = for i < ^, u}{t) = 1 for t > ^, (9) 



so uj'{x) = for X = 0,1. Then we set 



n-^ii 



-Xl\ 



which approximates the Green's function G associated with the domain 17* = (0, 1) x 



R ; in particular, it satisfies the boundary condition G 
the notations 



an. 



G 



9[d] ■■= g{d,X2,X3;$,), A 
we can rewrite the definition of G as 



6=0,1 

2axi/£ 



0. Using 



= [g[x^] -pg[-xi]] - [A 5[2-xi] - P >^'^ 9[2+xi]] ^(6)- 



(10) 



We now present a version of [5l Lemma 4.2], which gives certain upper bounds for 
g that involve a weight function A of type A^. 

Lemma 7. Let x £ [1 + e, 3] x M^. Then for the function g of ^ and the weight 
A := e^^^^i"^)/^ one has the following bounds 



UXg){x--)\U,n,<Ce, 
\\{Xd^,g){x--)h.n,<C, k = 1,2,3, 
\\{Xdlg){^;-)h.n,<Ce-\ k = 1,2,3, 



(11a) 
(lib) 
(11c) 



while for any ball B{x'; p) of radius p < g centred at x' € [0, |] x M? , one has 
ll(A5)(^;-)lli,i;n.nB(.';p) <Cmin{p/e, {p/ ef'^] e-^^'' . 



(lid) 



Proof. The bounds (11a), (lib) appear in [5l Lemma 4.2]. The estimate (11c) is 
shghtly sharper compared to the similar bound in [5] Lemma 4.2]; the latter is 
for the domain \B{x,p) and involves the logarithmic term ln(2 + e/p) as it is 
valid for xi G [1,3]. In the above Lemma [7] we make a stronger assumption that 
xi £ [1 + £,3], under which ^2* \ B{x,e) = $7* and the logarithmic term becomes 
In 3 so can be dropped. 

To obtain the final desired bound ( lid ) , note that xi>l + e implies that r > 1 
in so ([7a|) and ((Tb]) yield \ {\g\ + J^c/]) < C^-^a e"(«i-^) for k = 1,2,3. 
Here ii-?= -i\ii\+r) > -2l^i|so Ae"(«i-^) < qM^^i-^) / e-m ^ Note also that in 



B{x'; p) we have ^ 1 ~ g so ]^i] 



> + Consequently, A (J^J + 



|%5l) < C6-3Ae-i"/^ so ll(A5)(a;;-)lli,i;f7,nBK;p) < C{p/ef . Combining 

this with {p/ef < C{p/e)P el"'/^ for p = l,\ yields (|lld|. □ □ 

Lemma 8. Let il* := (0, 1) x M^, e G (0, cq] /or some sufficiently small constant cq, 
and < a < C T/ien the function G of (10) satisfies, for any x G [1/4,3/4] x M^, 
all i/ie bounds 0/ Lemma^ with g replaced by G. 
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Proof. Let D be any of the first- or second-order differential operators that appear 
in ([s]). Now for any fJ'^ C Jl*, the representation (10) yields 

\\DG{x; ■)\\i^ni = \\Dg[xi] - A"I)5[2-xi] +P^"^-Dff[2+xi]lli,ni 

-p||A+£'5[2+xi]l|i,ni. (12) 

For the first term that involves g[x-^], we use the corresponding lower bound from 
Lemma [5] so it remains to show that this bound will dominate the remaining three 
terms. For the second term we note that g[-xi] satisfies the upper bounds of type 
(jsl with O replaced by 0=,, [5]. Now, as xi > ^ implies that p < e~2"/^, the second 
term will be dominated by the first term if e is sufficiently small (i.e. if the constant 
Co is sufficiently small). 

To estimate the terms that involve X^g^2±xi] (12), we use the bounds (11) of 
Lemma [7] In particular, by (8d), (11c), 

l\\9lj^9{x;-)\\l-n,\B{x;p) - \\^^d'iig[2±xi]\\l,Q,\B{x;p) 

> e-^ [c ln(2 + e/p) - C] 

> ce~Mn(2 + e//)) 



for p < ci£ if ci is sufficiently small, so we get a version of (8d) for G. Similarly, by 
(8e), (11c), for k = 2,?, one gets 

\\\dl^9{x]-)\\l-n,\B{x;p) 



;fc5[2±xi]lll,Q,\B{cc;p) 

>e"^[c ln(2 + e/p) + c|lne| - C] 
> ce^'^{\n{2 + e/p) + |lne|) 



provided that e is sufficiently small. This yields a version of (Be) for G. Finally, 
(lid) implies that 

||(A±<7[2±x,])(a^;-)lli,i;f7.nB(x;p) < Cmin{/,/e,(p/e)V2}e-|°A 

< \c min{p/e, {p/ef/^} 

for any arbitrarily small c provided that e is sufficiently small. This observation 
yields a version of (pel) for G. □ □ 



Proof of Theorem^ for the domain $7* = (0,1) x . As, by Lemma [sj the 
approximation G satisfies the bounds that we need to prove for G, it suffices to 
estimate the function v = G — G, which satisfies the differential equation 

Clv{x;^) = -eA^v{x;$,) + 2ad^,v{x;$) = (l){x;^) for ^ G SI,, (13) 

and the boundary condition v{x; ^) \^^qq = 0. Here for the right-hand side (f), it was 
shown in [51, Lemma 5.1] that 



{x;-)\\i.n,<Ce-^^"/' 



(14) 
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for some constant C3. In view of (13), the function v can be represented using the 
Green's function G as 



G{s; ^) (f>{x; s) ds 



So applying to this representation with p = 0,1,2, k = 1, 2, 3, for any sub-domain 
il'^ C 17*, one gets 



mv{x;.)h,n',<i sup ||9|G(s;.)||i,nt 



< Ce-^W^- sup \\dlG{s;-)\\,^n,, 



sen 



where we also used (14). We shall use the above estimate with il^ := 0* for p = 0,1 
and := B(x, p) for p = 2. As the bounds ([s]) of Theorem [l] remain valid for 
the domain il* and also ||G(ic; •)||i,r2, < C* [5', 'G], one now concludes that 

Wdlvix; ■)h.n,\Bi.,p) < Ce^^^"/^ • e-\ln{2 + e/p) + \ lne\) 
< c'e"Mn(2 + e//9) for A; = 1,2,3, 

where c' is arbitrarily small provided that e is sufficiently small. Similarly 

<Ce-^^"/'-e-i/2 <c'e, 

which also implies 

\\v{x\-)\\i,i-n,nB{a,,p) < c min{p/e, (p/e)^/^}. 



(15a) 
(15b) 



(15c) 

As c' in the bounds (15) can be made arbitrarily small, combining them with 
Lemmajsjyields the desired bounds of Theorem|3]for the domain = (0, 1) xR^. □ 



4 Proof of Theorem 

(0,1)3 

The proof of Theorem [s] for the domain $7 = (0, 1)^ is very similar to the above 
proof for the domain Q,^: = (0, 1) x presented in Section [sj The only difference 
is that instead of the approximation G for the domain (0,1) x we now use the 
approximation G^ defined by 

Ga{x;^):=G{x;^) - coo{C^) G{x;Ci, -(2,(3) - u;i(6) G(^; 6, 2 - 6, 6), 
Gw{x;^) ■.= Ga{x;$)-iUo{^3)Gaix;^i,^2,-^3)-iOi{C3)Ga{x;Ci,^2,2-^3). 

Here u}o{t) ■= (^{1 — t) and uji{t) := uj{t) with uj defined in ([9]), so that 0Jk{k) = 1 
and ujk[l — k) = for k = 0, 1. This approximation was constructed employing the 
method of images; an inclusion of cut-off functions ensures that it vanishes on (90. 
All the properties of G given in Section [3] remain valid for this new approximation 
with Q.^ replaced by Q. provided that x G [j, |]^. We leave out the details and 
only note that the application of the method of images in the S,2- and ^3-directions 
is relatively straightforward as an inspection of ^ shows that in these directions, 
the fundamental solution g is symmetric and exponentially decaying away from the 
singular point. □ 



3 for the domain Q = 
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5 The Convection- React ion- Diffusion Case 

We now slightly generalize ([T]) by including a reaction term with a constant coeffi- 
cient /3 > 0: 

Cxu{x) = —e/S.xu{x) — 2adx{a{x) + j3u{x) = f{x) for x € Q (16a) 

u{x) = for X G dn. (16b) 

Now the fundamental solution g in satisfies, for each fixed a; G M'^ the following 
version of (5) with the adjoint operator 

^Igix; ^) = -eA^g{x; ^) + 2a %g(ic; + f3g{x; ^) = 6{x - ^) for ^ G Ml 
Again imitating a calculation of [21 [9] , one gets a version of (§: 

1 e"^i'[-il~'^''[="il 



9ix'^^) = j — 2 ^ ' where 7 := y^a^ + e/3. 

In view of 7 = a + + O (e^) , an inspection of the proof of Theorem [s] shows 
that the lower bounds Q remain valid for the convection-reaction-diffusion problem 



(16). 



6 Outlook for problems in n dimensions 

It was shown in [7j that the upper bounds of Theorem [l] remain valid for a two- 
dimensional variable-coefficient version of ([T]). Note that the fundamental solution 
g-^2 that solves ^ in M^, and also its derivatives involve the modified Bessel func- 
tions of second kind of order zero -K'o(') and order one Ki{-). This fundamental 
solution is given by 

ffM2(^;^) = ^e"«ii^o(af) 

with the notations x, ^, and r appropriately adapted. Using this explicit represen- 
tation, one can imitate the proof of Theorem [3] and get similar lower bounds for 
the two-dimensional case. A certain difficulty lies in having to deal with the Bessel 
functions, for which one can simply employ asymptotic expansions p!} I12|: 

Ko{z) = ( J)'^'e-^ (1 + {z-')) for \z\ » 1, 

Ki{z) = Ko{z)(l + ^ + (z-2)^ for |z| » 1, 

Ko{z) = -lnz + Oil), Ki{z) = ^ + 0{1) for |z| < 1. 

In this manner one gets the following result. 

Theorem 9. The Green's function associated with problem ([T]) in the unit-square 
domain Q := (0, 1)^, satisfies a two-dimensional version of Theorem^ 
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Finally, let us take a look at the problem ([T]) in the n-dimensional domain $7 := 
(0, 1)" of an arbitrary dimension n > 2. The corresponding fundamental solution 
5Rn is given by 

1 /a\"/2-i c 

= (2^)n/2gn-i \-^) e°^iK„/2-i(ar) 

with the modified Bessel functions i^n/2-i(") of second kind of (half-)integer order 
n/2 — 1, and the notations a;, ^, and r appropriately adapted. Using asymptotic 
expansions of these Bessel functions [HIE], one can again get a version of Theoremjsj 

References 

[1] M. Abramowitz and I. A. Stegun. Handbook of Mathematical Functions with 
Formulas, Graphs and Mathematical Tables. Applied Mathematics Series. Na- 
tional Bureau of Standards, Washingtion, D.C., 1964. 

[2] N.M. Chadha and N. Kopteva. Maximum norm a posteriori error estimate 
for a 3d singularly perturbed semilinear reaction-diffusion problem, doi: 
10.1007/sl0444-010-9163-2 (published online 1 June 2010), 2010. 

[3] K. Eriksson. An adaptive finite element method with efficient maximum norm 
error control for elliptic problems. Math. Models Methods Appl. Sci., 4:313-329, 
1994. 

[4] S. Franz and N. Kopteva. A posteriori error estimation for a convection- 
diffusion problem with characteristic layers, in preparation, 2011. 

[5] S. Franz and N. Kopteva. Full Analysis of Green's function estimates for a 
convection-diffusion problem with characteristic boundary layers in 3d. Tech- 
nical report, University of Limerick, 2011. 

[6] S. Franz and N. Kopteva. Green's Function Estimates for a Convection- 
Diffusion Problem in Three Dimensions, in preparation, 2011. 

[7] S. Franz and N. Kopteva. Green's function estimates for a singularly perturbed 
convection-diffusion problem, submitted, 2011. 

[8] J. Guzman, D. Leykekhman, J. Rossmann, and A. H. Schatz. Holder estimates 
for Green's functions on convex polyhedral domains and their applications to 
finite element methods. Numer. Math., 112:221-243, 2009. 

[9] R. B. Kellogg and S. Shih. Asymptotic analysis of a singular perturbation 
problem. SI AM Journal on Mathematical Analysis, 18(5): 1467-1510, 1987. 

[10] N. Kopteva. Maximum norm a posteriori error estimate for a 2d singularly 
perturbed reaction-diffusion problem. SIAM J. Numer. Anal., 46:1602-1618, 
2008. 

[11] R. H. Nochetto. Pointwise a posteriori error estimates for elliptic problems on 
highly graded meshes. Math. Comp., 64:1-22, 1995. 

[12] F.W.J. Olver, D.W. Lozier, R. F. Boisvert, and C.W. Clark. NIST Handbook 
of Mathematical Functions. Cambridge University Press, Cambridge, 2010. 



11 



